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Abstract 

The uncertainty in the equihbration in heavy ion colhsions due to the choice of the 
value of the couphng is examined. The results of the equilibration are indeed affected 
by this choice. In particular, a variation of as from 0.3 to 0.5 reduces the parton phase 
of the plasma by as much as 4.0 fm/c and increasing coupling causes a reduction in 
the generated entropy. As far as equilibration is concerned, larger coupling results in 
faster equilibration both chemically and kinetically but improvements only happen to 
the fermions. These are accompanied by more rapid cooling and therefore shortened 
lifetime. In the light of these results, to choose Og is almost equivalent to choosing 
the results. Also because of consistency, any fixed Os is incompatible with an evolving 
plasma. The best choice is therefore not to choose at all but let the system makes 
its own decision. We show a simple recipe how this can be done. With an evolving 
coupling, equilibration is accelerated with better results at the expense of the duration 
of the deconfined phase. 



1 Introduction 



In future heavy ion collision experiments at LHC and at RHIC, a major effort is not only to 
try to produce deconfined matter, given some of which might have already been created at 
CERN SPS at present energies, but also to show beyond any doubt that during the collisions 
that deconfined matter really exists for however brief moment. At AGS and SPS, this latter 
task would be difficult because their energies might not be sufficiently high or the system 
might not be large enough to reduce boundary effects on the produced quark-gluon plasma. 

In the very violent collisions, many particles will be produced, some of which can escape 
from the system such as electromagnetic probes, others can only reveal themselves after the 
final break up of the system such as probes using strange and charm hadrons. Nevertheless, 
they will all be affected by the evolution of the system. Particle production are not exactly 
the same during the initial, pre-equilibrium and equilibrium phases, therefore the durations 
of these periods also play a role in modifying the yield of the final produced particles which 
will eventually fall into the detectors. It is therefore necessary to study the evolution of the 
system or equivalently to study the equilibration of the parton plasma. There have already 
been quite a few works on the equilibration in the parton phase of the plasma, for example 
chemical equilibration ^\ thermalization '^^ and full equilibration Our purpose here is 
not exactly to repeat these investigations but to concern ourselves on one of the uncertainties 
which can alter the results of the equilibration. The first and most obvious is caused by the 
uncertainty in the initial inputs. These can be studied but we will not do this at present. 
The second uncertainty could have been the infrared screening parameter used in the usual 
perturbative QCD. However, in a multi-particle system, QCD will generate Debye, quark and 
gluon medium masses ^\ These will effectively screen off the infrared divergences. So this 
parameter and therefore the uncertainty needs not be present. The third one is the strong 
coupling constant itself. It is the only remaining free parameter once that of infrared screening 
is no longer present. We would like to find out how the results will change with the value 
of as- There is however a second reason concerning consistency for which we would like to 
look into the uncertainty arising from the choice of a^. The common choice for these kinds 
of studies is as = 0.3 which corresponds to an average momentum transfer of Q ~ 2.0 GeV 
with Aqcd = 200 MeV. This choice is reasonable early on after the initial collisions. We 
show in Eq. (|l|), the evolution of the average energy for quarks and gluons in our previous 
investigation into the equilibration at LHC and at RHIC. As can be seen, the average parton 
energy drops by at least 1.0 GeV over the entire evolution. Since the average momentum 
transfer has to be related to the average parton energy, and so as = 0.3 cannot be a good 
value for the entire duration of the evolution. In fact, any choice of fixed will equally not 
be good enough. We will show a recipe for obtaining a varying coupling in the next section 
to overcome this problem of consistency and will show the corresponding results as well as 
those with other fixed values of the as- Together they reveal the effects of the coupling on 
the equilibration. 
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Figure 1: The average gluon (solid) and quark (dashed) energies change a lot during the 
evolution of the parton plasma both a LHC and at RHIC. Therefore the average momentum 
transfer should also vary considerably and hence the coupling as well. 

2 Evolution of a parton plasma with various couplings 

To study and to make the effect of as manifest, we choose some large values of as = 0.5, 
0.8 to do the evolution in addition to the previous as = 0.3 case. Now also to overcome the 
problem of consistency, we introduce, as a solution, a varying coupling which evolves with the 
system in the following way. Since two incoming partons each carrying the average parton 
energy {E) = {ttot) I {^tot) can exchange a maximum = A{E)'^. So combining the {E{t)) at 
any moment r and the 1-loop running coupling formula, we have a coupling that is entirely 
determined by the system. With this latter approach, there is no remaining free parameter. 
This is in fact, with hindsight, a better choice for the coupling as we will see presently. We 
denote this coupling by a^ from now on. 

With these choices of the coupling, we show the effects of the coupling in Fig. ^ and 
Fig. 1^. Fig. § shows the effects on chemical equilibration in terms of the fugacities as well as 
the parton estimated temperatures. Fig. ^ shows the pressure and energy density to pressure 
ratios from which we can deduce information on kinetic equilibration. Let us first look at 
parton chemical equilibration. We used the same fixed initial conditions as before so as 
to concentrate only on the coupling. From Fig. ^, we see that with increasing a^, chemical 
equilibration is definitely faster both for gluons and for quarks, however, only in the case 
of the fermion, do they show any improvements. For gluons, the end degree of chemical 
equilibration does not change very much with the coupling. 

For kinetic equilibration, we use the pressure to pressure and energy to pressure ratios as a 
check of the isotropy of the parton momentum distribution. Because of our thermalized initial 
conditions, these ratios should start at 1.0 or with an isotropic distribution. This isotropy is 
lost subsequently, as seen in Fig. ^, where all the curves shift downward away from 1.0. This 
is due to the disruptive effect of the longitudinal expansion. As the net interaction responds 
by increasing its rate, the expansion effect is later overcome and isotropy is progressively 
being recovered. This is when all the curves in Fig. |] rise again. Larger couplings lead to 
faster equilibration but as in the case of chemical equilibration, only quarks and antiquarks 




Figure 2: With increasing coupling, chemical equilibration is faster but the cooling is also 
more rapid. As seen here, the fastest increase in the parton fugacities Ig and Iq are the 
curves (dashed) with = 0.8, the next are those produced with as = 0.5 (dotted) and 0.3 
(solid). The curves of (long dashed) shift across the constant as "contours" and equihbrate 
definitely better than the solid lines. The shifts of the descending T curves with as are sizable. 
As a result, the lifetime of the parton phase of the plasma is controlled by the value of the 
coupling. 

show obvious improvements. Equilibration is clearly faster for all partons in Fig. |^ and Fig. ^ 
However, this is achieved at the expense of similarly faster cooling. The temperature estimates 
all decrease more rapidly than before in Fig. |^. Not only that but the time at which they 
reach the assumed phase transition temperature at Tc ~ 200 MeV changes by as much as 4.0 
fm/c when the varies from 0.3 to 0.5 at LHC So the lifetime of the parton phase of the 
plasma is very sensitive to the value of as- Similarly sensitive to the coupling is the generated 
entropy It is clear therefore, as we have already mentioned, to have to choose a value for 
as is not the best choice. With the consistent a^, we see in Fig. ^ and Fig. ^ these curves 
start off near as = 0.3 case but shift away across the constant as "contours". They therefore 
achieve faster and better equilibration for quark and antiquark than the as = 0.3 case. They 
equilibrate, in this case, in an accelerated fashion. This is a non-abelian effect not found in 
the equilibration of ordinary electromagnetic plasma or in other many-body system. 

The non-abelian effect mentioned in the previous paragraph can be shown more clearly 
via the collision time 9 defined previously through the relaxation time approximation ^^^^ . In 
Fig. we have plotted the evolution of these for the various values of the coupling. From 
the fixed coupling 6''s, their behaviours are similar, i.e. a fast rapid initial decrease and 
then a slow rise until the end. These can be explained by the short expansion and the net 
interaction dominated phases. It would be helpful to identify the inverse collision time as 
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Figure 3: These ratios show the isotropy of the parton momentum distribution and therefore 
kinetic equihbration with r. The top (bottom) set of four curves are for energy (pressure) to 
pressure ratios. Faster thermahzation is seen with increasing couphng. The assignment of the 
couphng to the curves are = 0.3 (sohd), 0.5 (dotted), 0.8 (dashed) and (long dashed). 
Improvements are, however, reserved only for the fermions. 
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Figure 4: Evolution of the collision time reflects that of the net interaction rate. The values 
of the coupling are assigned to the curves in the same way as before. The curves produced 
with fixed couplings have similar time-dependent behaviours. The case of is, however, very 
distinct. It shows the non-abelian effect of QCD accelerates the equihbration. This is unique 
to a QCD parton plasma. 



the net interaction rate, then the first phase is simply the response of the system to being 
driven out of equihbrium. As a consequence, the net interaction rate has to increase or 9 
has to drop. As equihbrium is approached, the net rate has to slow down so 6 has again to 
increase but slowly. In the case of a^, the slow rise of the net interaction dominated phase is 
replaced by a continued slow decrease. So the net rate continues to increase albeit at a much 
slower rate as the initial drop. We have already seen in Fig. ^ and Fig. |^ that equilibration 
is present in all cases, so what is happening to the a1 case is, the increase in the interaction 
strength compensates for the near equilibrium slowing down of the net rate. The result is the 
non-abelian accelerated equilibration that we have already seen. This is of course unique to 
a QCD plasma. 

In summary, we have studied the dependence of the equilibration and introduced a 
simple recipe to solve the consistency problem raised in the introduction. It is found that 
larger means faster equilibration for all partons but improvements are only for the fermions. 
Lifetime of the parton phase as well as the entropy are sensitive to the value of as therefore 
the consistent is a better choice which gives rise to accelerated equilibration unique to a 
QCD plasma. 
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